Statistical mechanics in the context of special relativity II 
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The special relativity laws emerge as one-parameter (light speed) generalizations of the corre- 
sponding laws of classical physics. These generalizations, imposed by the Lorentz transformations, 
affect both the definition of the various physical observables (e.g. momentum, energy, etc), as well 
as the mathematical apparatus of the theory. Here, following the general lines of [Phys. Rev. E 
66, 056125 (2002)], we show that the Lorentz transformations impose also a proper one-parameter 
generalization of the classical Boltzmann-Gibbs-Shannon entropy. The obtained relativistic entropy 
permits to construct a coherent and selfconsistent relativistic statistical theory, preserving the main 
features of the ordinary statistical theory, which recovers in the classical limit. The predicted dis- 
tribution function is a one-parameter continuous deformation of the classical Maxwell-Boltzmann 
distribution and has a simple analytic form, showing power law tails in accordance with the exper- 
imental evidence. Furthermore the new statistical mechanics can be obtained as stationary case of 
a generalized kinetic theory governed by an evolution equation obeying the H-theorem and repro- 
ducing the Boltzmann equation of the ordinary kinetics in the classical limit. 

PACS numbers: 03.30.+p, 05.20.-y 



I. INTRODUCTION 

In high energy physics, various approaches of statistical 
mechanics and kinetic theory have been adopted to study 
cosmological models 0, 0, multiparticle production in 
hadron-hadron collisions [3(, possible Lorentz invariance 
violation in extension of standard- model , neutron star 
matter black hole models [EH 0,0 etc. Common to 
all such approaches is the explicit or implicit assumption 
of a very specific form of the entropy, i.e. the celebrated 
Boltzmann-Gibbs-Shannon entropy which, according to 
Jaynes |lOj maximum entropy, principle conducts to the 
Maxwell-Boltzmann exponential distribution. 

On the other hand since long time it is known that 
there exist open problems in high energy physics which 
cannot be solved within the ordinary Maxwell-Boltzmann 
statistics e.g the cosmic rays spectrum [H[ll El and the 
black hole entropy Bekenstein-Hawking area law [T3.ll5| . 

In statistical physics, beside the ubiquitous exponen- 
tial distribution cxp(— x), also non-exponential distribu- 
tions have been considered. The most famous ones ap- 
pear within the quantum statistical mechanics, namely 
the Fermi-Dirac and the Bose-Einstein distributions, 
which are forced by the exclusion-inclusion principle. In 
the last decades a particular attention has been addressed 
to the statistical distributions presenting power law tails 
as x~ a . The latter distributions have been observed ex- 
perimentally in man y fi elds of physics, including cosmic 
rays physics 0,0, El; plasma physics 0] and multi- 
particle production processes |l7j . 

In order to propose theories based on non-exponential 
statistical distributions, one can follow two different 
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ways. The first is the entropic one consisting of a proper 
generalization of the system entropy. After fixing the en- 
tropy, the Jaynes maximum entropy principle univocally 
determines the form of the distribution function. The 
problem then is essentially reduced to the choice of the 
right entropic form. 

The second way is the kinetic one which is very fa- 
miliar in non equilibrium statistical mechanics. It is 
well known that the classical statistical mechanics can 
be obtained as particular case of the ordinary kinetic 
theory for t — > oo, within the Boltzmann or the Fokker- 
Planck picture. Generalized kinetic theories necessarily 
conduct to generalized statistical mechanics producing 
non-exponential distributions. 

Both the entropic as well as the kinetic path, conduct 
to the conclusion that it is possible to develop, beside the 
classical statistical mechanics, other coherent and self- 
consistent statistical theories [ll llH Efl Ell 151 151 Ety . 

It is remarkable that the Maxwell-Boltzmann distri- 
bution does not emerge within the statistical mechan- 
ics. It is more proper to say that the classical statistical 
mechanics is built starting from the Maxwell-Boltzmann 
distribution. It is commonly accepted that this distri- 
bution emerges within the Newton mechanics. Indeed, 
numerical simulation of classical molecular dynamics un- 
equivocally conducts to the exponential distribution. At 
this point, the question if the Maxwell-Boltzmann dis- 
tribution obtains also in the case when the microscopic 
dynamics is governed by the special relativity laws, nat- 
urally arises. 

In the relativistic statistical theory the collision invari- 
ants are imposed by the laws of special relativity. Hence 
the first difference between the classical and the rela- 
tivistic distribution functions regards the different forms 
of the collision invariants, in the argument of the dis- 
tribution function. A second difference could be origi- 
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nated by the form of the entropy of the relativistic many 
body system, which could be different with respect to the 
Boltzmann-Gibbs-Shannon entropy. 

In the ordinary relativistic statistical mechanics it is 
accepted that the entropy has the same form of the 
classical statistical mechanics. Then the Jaynes princi- 
ple yields an exponential distribution which reproduces 
exactly the Maxwell-Boltzmann distribution in the rest 
frame. The most famous and important relativistic par- 
ticle system is undoubtedly the cosmic rays 0, 0, ^1 . 
The cosmic rays spectrum has a very large extension 
(13 decades in energy and 33 decades in particle flux) 
and presents a power law asymptotic behavior. Unfor- 
tunately the Maxwell-Boltzmann distribution fails to ex- 
plain this experimental spectrum. The same power law 
asymptotic behavior of the distribution function was ob- 
served also in other relativistic systems ^| . Hence 
for relativistic particle systems, the experimental evi- 
dence suggests a non-exponential distribution function 
with power law tails. This distribution, can be origi- 
nated exclusively by an entropy manifestly different from 
the Boltzmann-Gibbs-Shannon one. 

We recall that all the classical physical observables 
(momentum, energy, etc) when considered within the 
special relativity are properly generalized. Any rela- 
tivistic formula can be viewed as a one-parameter (light 
speed) generalization or deformation of the correspond- 
ing classical formula. Consequently it is natural to as- 
sume that also the entropy of a relativistic system could 
be a one-parameter generalization of the classical en- 
tropy. In this case the logarithm appearing in the expres- 
sion of Boltzmann-Gibbs-Shannon entropy must be re- 
placed by another one parameter generalized logarithm. 
Then also the ordinary exponential function appearing 
in the Maxwell-Boltzmann distribution must be replaced 
by an one-parameter generalized exponential. These two 
generalized functions have been proposed heuristically in 
ref. [Tsj and are given by 

ex P{K} (t) = (Vl + ^ 2 + Ktj , (1.1) 

\ny(t) = ^^- (1-2) 

Note that in the classical limit k — ► one easily re- 
covers the ordinary exponential and logarithm namely 
ex P {0} (t) = exp(t), ln {0} (t) = ln(t). Successively in ref. 

fli^l it has been shown that the statistical theory, devel- 
oped starting from these generalized functions, emerges 
within the framework of the special relativity. The en- 
tropy for a relativistic many body system assumes the 
form 

S (f) = -J d "P f ln {M(/) = - < h V>(/)> . (1-3) 

while the relevant statistical distribution in the rest frame 
is given by 



W being the relativistic kinetic energy, while a and A are 
two constants depending on k. These two latter relation- 
ships are one-parameter generalizations of the Bolzmann- 
Gibbs-Shannon entropy and of the Maxwell-Boltzmann 
distribution respectively, which are recovered in the clas- 
sical limit k — ► 0. 

Main goal of the present contribution is to show that 
the functions exp (i) and ln {K} (f), leading to the dis- 
tribution function (|1.4|) and to the entropy (|1.3|) . emerge 
naturally and unequivocally within the special relativity 
theory. These functions replace the ordinary exponential 
and logarithm, relegating them at the native status of 
classical functions. It is remarkable that the relativistic 
statistical theory, based on the entropy (|1.3fl . predicts a 
power law asymptotic behaviour to the spectrum of rela- 
tivistic particles according to the experimental evidence. 

The paper is organized as follows: 

In Sect. II we recall some concepts of relativistic dy- 
namics focusing our attention on the deformations in the 
mathematical formalism enforced by the finite value of 
the light speed. 

In Sect. Ill we show that the form of exp {K} (t) can be 
obtained within the special relativity by using elemen- 
tary algebraic calculations, starting from the dispersion 
relation of free relativistic particles. 

In Sect. IV we reobtain the exponential exp {K} (i) start- 
ing from the concept of Lorentz invariant integration. 

In Sect. V starting from the three-dimension K-sum 
enforced by the Lorentz transformations we introduce the 
concept of generalized differential. 

In Sect. VI we derive the forms of ln {K} (i) and exp {K} (t) 
within the special relativity, starting from the additivity 
law of relativistic momenta. 

In Sect. VII we recall some concepts related to maxi- 
mum entropy principle, in the context of a class of gen- 
eralized statistical theories. 

In Sect. VIII we consider the explicit form of the dis- 
tribution function within the special relativity. 

In Sect. IX we consider the main features of the sta- 
tistical mechanics and of the kinetic theory based on the 
entropy (|1.3|) . 

In Sect. X we consider some applications of the present 
theory to specific physical problems and compare the re- 
sults with the ones obtained within the classical statisti- 
cal mechanics. 

Finally in Sect. XI some concluding remarks are re- 
ported, while the main mathematical properties of the 
functions exp (t) and ln, K ,(£) are collected in Appendix. 



II. DEFORMED SUMS IN SPECIAL 
RELATIVITY 

In this section we introduce the special relativity in an 
alternative way, starting from the well known formulas 
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of hyperbolic trigonometry (pag. 58 of ref. |2E 



arcsinh x — arccosh yl + x 2 = arctanh 



(2.1) 



We pose x = nq, being q the dimensionless variable cor- 
responding to the momentum and n a dimensionless de- 
formation parameter. After introducing the functions 



u(q) 



y/1 + K 2 q 2 



(2.2) 



(2.3) 



corresponding to the dimensionless velocity and kinetic 
energy respectively we can write Eq. H2.1(l in the form 

arcsinh(Kq) = arccosh(l + k 2 W) = arctanh(Ku) , (2.4) 

which in the classical limit k — > reproduces the dis- 
persion relation q = \J2~W = u of Newton mechanics. 
Finally, after introducing the dimensionless total energy 

£(W)=W + 4 . (2-5) 

one can write Eq. I|2.1|l also in the form 

arcsinh(Kq) = arccosh(K 2 £) = arctanh(Kw) . (2.6) 

It is immediate to verify the validity of the following for- 
mulas 



q(u) 

q(W) = \/2W + k 2 W 2 



\/T— K 2 U 2 



q(£) = 

u(W) 

u{£) = 
W(u) 



V2W + k 2 W 2 

1 + K 2 W 

^Jn 2 £ 2 - 1/n 2 

k 2 £ 
1 



K 2 \/l — K 2 U 2 
1 



W(£) = £-- , 



£{q) = -^l + n 2 q 2 , 

£{u) = 



K 2 \J\ — K 2 U 2 



(2.7) 
(2.8) 
(2.9) 

(2.10) 

(2.11) 
(2.12) 

(2.13) 
(2.14) 
(2.15) 



At this point one can define the physical variables ve- 
locity v, momentum p and total energy E through 



v p 
u mq 



kc = u,„ 



(2.16) 



and the kinetic energy as W = E — mc 2 , being c the light 
speed. We impose that lim^oo. v* < oo, in order 



to preserve the validity of the definitions of the physi- 
cal variable, also in the classical limit. After inserting 
the physical variables in the above obtained formulas, 
one recover all the formulas of the one particle dynamics 
within the special relativity. For instance, Eq.s l|2.7l) and 
(|2.14|) transform 



P 



E = y m 2 c 4 + p 2 c 2 



(2.17) 
(2.18) 



and so on. 

Let us consider in the one-dimension frame S two iden- 
tical particles of rest mass m. We suppose that the first 
particle moves towards right with velocity u\ while the 
second particle moves towards left with velocity u-i. The 
momenta of the two particles are given by q\ = q (ui) and 
q2 = q (1*2) respectively, where q (u) is defined through 
Eq. (|2.7I) . The total energies of the two particles are 
given by £\ — £ (u\) and £2 — £ (1*2) respectively, with 
£ (u) given by Eq. (|2.15|l . Analogously the kinetic ener- 
gies of the two particles are given by Wi = W (iti) and 
W 2 = W (112) respectively, with W (u) defined by Eq. 
(12321 . 

We consider now the same particles in a new frame 
S ' which moves at constant speed U2 towards left with 
respect to the frame 5. In this new frame, which is the 
rest frame of the second particle, the two particles have 

K 

velocities given by u ^ = Ui © 112 and u ' 2 = respectively, 
being 



U\ ® 1*2 



U\ + U2 
1 + K 2 U\U2 



(2.19) 



the well known relativistic additivity law for the dimen- 
sionless velocities. In the same frame iS ' the particle 
momenta are given by q[ = q(u'-y) and q' 2 = 0. Anal- 
ogously in S', the particle kinetic energies are given by 
Wi = W(u[) and W 2 = 0, while the total particle 
energies results to be £ [ — £ (u^) and £' 2 — l/n 2 . 

A very interesting result follows from the relation £' 2 = 
1/k 2 , regarding the physical meaning of the parameter 
k. It is evident that the quantity 1/k 2 represents the 
dimensionless rest energy of the relativistic particle. 

Let us ask now the following questions: is it possible 
to obtain the value of the momentum q [ (or of the total 
energy £ [ , or of the kinetic energy VV [ ) starting directly 
from the values of the momenta q\ and q2 (or of the total 
energies £\ and £ 2 , or of the kinetic energies W 1 and 
W 2) in the frame 5? The answers to the above questions 
are affirmative as we will see in the following. 

First we consider the case of the relativistic momentum 

K 

q[. Clearly we can write q[ = q(u[) = q{u\®U2)- Then 
following ref. we have 



q(ui © u 2 ) 



u\ © U2 



1 - K 2 (U1 © U2) 2 
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Ml © U 2 



\/l - K 2 (-rr 1 P :2 -) 2 

1 + K 2 U\U 2 



(til © u 2 



^(l-/i 2 U 2 )(l-K 2 U 2 ) 
til + «2 



^(l-K 2 li 2 )(l-K 2 W 2 ) 
"1 /-, 



9 2 



"2 / « 2 ^1 



+ 

= g^/l + K 2 g§+<fc^/l + /« 2 9? . (2.20) 
Hence we can write 

q (mi © u 2 ) = © ?M , (2.21) 

being 



92 = (?ia/1 + k 2 (z| + g2\/l + K 2 Q , i 



(2.22) 



the K-sum of relativistic momenta introduced in ref. |18| . 
In words, the relativistic momentum q [ of the first par- 
ticle, in the rest frame of the second particle S' , is the 
K-deformed sum of the momenta q± and q 2 of the two 
particles in the frame S. The K-sum of the relativistic 
momenta and the relativistic sum of the velocities are 
intimately related. They reduce to the standard sum as 
the velocity c approaches to infinity (or equivalently the 
parameter n approaches to zero). The deformations in 
the above sums, in both cases, are relativistic effects and 
are originated from the fact that c has a finite value. 

We proceed now to calculate the total energy £ [ of the 
first particle in the frame S ' starting directly from the 
values, in the frame 5, of the total energies £ \ and £ 2 . 
Preliminarily we recall the following mathematical prop- 
erties of the inverse hyperbolic trigonometric functions 
(pag. 58 of ref. ||) 

arcsinh X\ ± arcsinh x 2 



= arcsinh \x\\J\+x\ ± x 2 \j l + x\ ), (2.23) 
arccosh y\ ± arccoshj/2 



= arccosh[ 2/ i 2/ 2±V(y 2 -l)(y 2 2 -l) ) » (2-24) 



arctanh z\ ± arctanh z 2 

z\ ± z 2 



— arctanh ■ 



1 + Z\Z 2 



(2.25) 



After posing x = nq, y — k 2 £, z — ku and taking into 
account Eq. H2.6|) , we are able to define the k-suul of the 
relativistic total energies as follows 



£ 1 ®£ 2 = n 2 £ 1 £ 2 - 



^^£ 2 -l)(^£ 2 -l) .(2.26) 



After having noticed that 

k 1 1 K 1 1 

£{u x ) © £{u 2 ) = — © — 



1 



1 



y/(l-K 2 ui)(l-K?ui) 



1 — K 2 u\ 



- 1 



1 — K 2 U 2 



- 1 



1 

1 




^ V(l - 


- K 2 U 2 )(1 — K 2 U 2 ) 
UlU 2 




2 U 2 ){l-K 2 U 2 ) 


1 


1 + K 2 U\U 2 




- K 2 U 2 )(1 — K 2 U 2 ) 


1 


1 + K 2 U\U 2 


« 2 


f K 2 UlU 2 ) 2 — K 2 (ui + U 2 ) 2 


1 


1 


K 2 1 

v 1 - 


r 2 ( u 1 +u 2 \ 2 

^ 1 + K 2 «ltl 2 J 


i 


1 



( ui © u 2 ) 



(2.27) 



we can conclude that the total energy £ [ = £{u\ © 
M2) of the first particle, in the rest frame of the second 
particle S', is the K-sum, defined through Eq. (|2.2tj[) . of 
the total energies £\ and £ 2 of the particles in the frame 
S: 

£ (ui © u 2 ) = £[ux) © £ (tta) . (2.28) 

The K-sum for the relativistic kinetic energies can be 
defined starting from 



m © W 2 = £ 1 © £ 2 - \ 



(2.29) 



and assumes the form 



Wi © W 2 = Wi+W 2 +k 2 WiW 2 

+ y / WiW 2 (2+«: 2 >Vi)(2 + K 2 >V2) ■ (2.30) 

After taking into account Eqs. (|2.12() . (|2.19|) and (|2.30|) 

one immediately obtains 



W fui © u 2 ) = W(ui) © W(u 2 ) 



(2.31) 



In words the K-sum, defined through Eq. (|2~5U|l . of the 
particle kinetic energies Wi and W 2 in the frame 5, gives 
the kinetic energy of one of the two particles in the rest 
frame S ' of the other particle. 
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Besides the relationships (|2.21|l . (|2.28|) and (|2.31|l . link- 
ing the K-sums defined through Eqs. lj2~T$j) . jg^ , (l2~2lj|) 
and H2.30[l . one can easily obtain: 



© g(u 2 ) = g ( tti 8 w 2 J , 
© £{u 2 ) =^Uil u 2 ) 
VV(ui) © W(u 2 ) = W Ui © u 2 

K / K 

u(qi) © u(q 2 ) = u ( qi © g 2 
£(<?i) ®^(g 2 ) =£(g<i ©g 2 
V%i) © W(g 2 ) = Wke g 2 

®q{£ 2 ) =q(£i ®£% 
u{E l )®u{E 2 ) = u(s 1 ®£ 2 
W{£i) © W(£ 2 ) = W (fi © £ 2 
q(Wl) ®q(W 2 ) =q(Wi © W 2 
© u(W 2 ) = u (>Vi © W 2 
£{Wi) © £ (W 2 ) = £ (Wi © W 2 



(2.32) 
(2.33) 
(2.34) 
(2.35) 
(2.36) 
(2.37) 
(2.38) 
(2.39) 
(2.40) 
(2.41) 
(2.42) 
(2.43) 



Let us make some remarks on the meaning of the k- 
sums. We note that in the classical limit the K-sum for 
the relativistic kinetic energies given by Eq. (|2.30(l re- 
duces to the following expression of Galileian relativity 



Wi © W 2 = Wi + W 2 + 2V>ViW 2 



(2.44) 



Then a deformed sum appears also in classical physics 
when we change the frame of observation of the particle 
system. Clearly this happens only for the kinetic energy. 
For the classical momenta and velocities hold the ordi- 
nary sum. 

We emphasize that the four K-sums given by Eqs. 
(|2~Tg|) . ^TF^i . (|2~2T)|) and (j2~3Uj) . emerge only when we 
change the frame of observation of the particles from S 
to S' . More precisely these sums give the value of q[ (or 
u'x, £'x, VV'i) of the first particle in the rest frame S 1 of 
the second particle, starting directly from the value of q\ 
(or u\,£\, Wi) and q 2 (or u 2 ,£ 2 , W 2 ) of the two particles 
in the old frame S. 



III. MINKOWSKI FOUR- VECTORS AND 
k-EXPONENTIAL 

From the general discussion in the previous section ap- 
pears clear that any relativistic formula can be viewed as 
a one-parameter deformation of the corresponding clas- 
sical formula. For instance the kinetic energy given by 



the relativistic dispersion relation 

W = \J m 2 c i + p 2 c 2 — mc 2 



(3.1) 



can be viewed as the one-parameter relativistic deforma- 
tion of the classical parabolic dispersion relation W c i = 
p 2 /2m. In the present section we will show that the rel- 
ativistic dispersion relation imposes a one-parameter de- 
formation to the ordinary exponential function, which 
presents now as a mathematical tool of classical physics. 
The new deformed exponential, as we will see in the fol- 
lowing sections, replaces the ordinary exponential in the 
relativistic physics. 

We adopt the metric — diag(l, —1, —1, —1) in the 
Minkowski space and recall that the length of any four- 
vector is Lorentz invariant. In particular for the four- 
momentum p^ = (E/c, p), p M = {E/c, — p) one obtains 
the relativistic dispersion relation 



„.2 2 

m c 



(3.2) 



which in terms of the total energy assumes the form 



E = y / m 2 c i + p 2 c 2 



We can write the latter equation as follows 

2 



E 



-m =i 



and then 



(* 


--) 






\ mc 


mc J 


\ mc 


mc I 



(3.3) 



(3.4) 



(3.5) 



At this point we can eliminate the variable E by employ- 
ing once again Eq. I|3.3|l obtaining 




(P_\ _ P 



\mc/ 



DIC 




p y v_ 

mc) mc 



1 .(3.6) 



The latter relationship, after introducing the dimension- 
less momentum q, defined through Eq. 12.16|) , can be writ- 
ten as 

(y/l + K 2 q 2 - Kq^j Ul + tfq* + Kq) = 1 , (3.7) 

or equivalently 

(yi + K 2 q 2 - Kq) 1 * (yi + K 2 q 2 + Kq) lA = 1 . (3.8) 

We remark that Eq. (|3.8(l follows directly from the dis- 
persion relation. Interestingly in the classical limit n — > 
Eq. I|3.8|) reduces to exp(— q) exp(q) = 1 while the dis- 
persion relation becomes the classical one W = p 2 /2m. 
In this way we obtain a direct link between the disper- 
sion relation of free classical particles and the ordinary 
exponential function. 
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In the light of the above result, we reconsider Eq. (|3.8() 
which, after introducing the function 



1/re 



ex P {K }(?) = (yi + KV + K-q 
writes in the form 

ex P {K} (-<7) exp {K} (g) = 1 . 



(3.9) 



(3.10) 



The function exp {K} (q), reproducing the exponential 
function, exp(g) = ex P {0} (?) in the classical limit re — > 0, 
represents a one-parameter relativistic generalization of 
the ordinary exponential. 



IV. LORENTZ INVARIANT INTEGRATION 
AND k-EXPONENTIAL 

In this section we present a second path which un- 
ambiguously conducts to the exp (q) function. Let us 

{K} 

consider the following integral in the three dimension mo- 
mentum space, within the classical physics framework 



d 3 qF 



(4.1) 



being F an arbitrary quantity. Whenever F depends 
only on q — |q|, the above integral can be reduced to the 
following one dimension integral 



I cl = / dqAnq z F . 

'o 



(4.2) 



After introducing the physical variable p = mv%° q the 
integral (|4.1|l transforms 



Id = 



d 3 p 



(mi> J° ) ■ 



F 



(4.3) 



In the framework of a relativistic theory it is well 
known that the integral (|4.3|) must be replaced by a 
Lorentz invariant integral I c i — ► I re i, being 

I rel =A fd 4 p 9(p ) 5(p% - m 2 c 2 ) F , (4.4) 



and A a constant. Now we can introduce the variable 
q = p /mu* and take into account that p^ — (E/c, p) 
and E = y m 2 c 4 + p 2 c 2 . After properly choosing the 
value of A, the four-dimension integral l|4.4|l can be easily 
reduced to the following three dimension integral 



d 3 q 



K 2 q 2 



F 



(4.5) 



We remark that in l|4.4|) the integration element d A p is a 
scalar because the Jacobian of the Lorentz transforma- 
tion is equal to unity. Then I re i transforms as F. For this 
reason in 14.5( 1 the integration element d 3 q/ yl + n 2 q 2 is 



a scalar. Finally we can reduce l|4.5(l in the following one 
dimension integral 



Ifpl — 



dq 



n 2 q 2 



Airq F 



(4.6) 



It is important to note that in classical limit re — > 0, Eqs. 
(|4.5|) and (|4.t)|) reproduce the corresponding classical ones 
given by (|4.1() and (|4.2() respectively. 

We focus now our attention to the classical and rel- 
ativistic expression of the one dimension integrals given 
by (|4.2I) and 1)4. respectively. One immediately observes 
that the relativistic integral is obtained directly from the 
classical one, by making the substitution dq — > d {K} q, 
being 



d {K} q 



dq 



v/l + n 2 q 2 



(4.7) 



the ^-differential. The relativistic one dimension integral 
can be written in the form 



poo 

Irei = / d {Ki q 4nq 2 F 
Jo 



(4.8) 



Clearly the K-integral J dq {K} is originated from the fact 
that the light speed has a finite value. 

As a working example, involving explicitly the re- 
integration, we consider the definition of the kinetic en- 
ergy. In classical mechanics the kinetic energy, using di- 
mensionless variables, is defined as 



W c i(q)= [ dqq 
Jo 



(4.9) 



It is straightforward to verify that the extension of the 
above definition in special relativity is given by 



[ d {K} q q 
Jo 



(4.10) 



Then, the replacement of the ordinary integration with 
by the K-integration, in the classical definition l|4.9|l . per- 
mits us to recover the relativistic expression l|2.3|l of the 
kinetic energy. 

Let us consider now the re-derivative related to the 
above defined re-integration 



d {li} q dq 



(4.11) 



Eq. H4.1Ufl linking W and q can be written now also in 
the following differential form 



d {K} q 



V%) = q 



(4.12) 



which, after integration with the condition W(0) = 0, 
yields the relativistic expression of the kinetic energy 
(|2.3|) . In the re — > limit the above differential equa- 
tion reduces to the classical one (d/dq)W(q) = q. 
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Taking into account that the ordinary exponential 
results to be eigenfunction of the ordinary derivative, 
namely [d/dq) exp(g) = exp(q), the question to deter- 
mine the eigenfunction of the K-derivative, naturally 
arises. After some simple calculations one obtains 



ex P {M («) = ex P {M (?) > 



(4.13) 



so that the K-exponential arises as eigenfunction of the 
K-derivative. 



V. RELATIVISTIC SUMS AND 
k-DIFERENTIAL 

In this section we explain the physical origin of 
the mechanism enforcing the deformation in the k- 
differential. We consider two identical relativistic par- 
ticles with rest mass m and velocities vi and V2 in the 
three-dimension frame S. The modulus of the relative 
velocity V of the particles depends on Vi and V2 and is 
given by (page 20 of ref. [26|) 



V(vi, v 2 ) 
being 



(vi e v 2 y 



vi e v 2 = 



1 



Vi X v 2 
1 - V!V 2 /c 2 



Vi - V 2 



1 - viv 2 /c 2 



(5.1) 



(5.2) 



We perform now the calculation of the modulus of rel- 
ative momentum related to the modulus of the relative 
velocity of the two particles according to the relativistic 
formula 



P(V) = 



iV 



y/1 - V 2 /c 2 



(5.3) 



Clearly results that P(V) 
into account that 



-P(vi,v 2 ) and after taking 



p/m 



v/T+ p 2 /m 2 c 2 



(5.4) 



we can conclude that P(V) = f(pi,Pi)- At this point 
we introduce the dimensionless variables q = p/rav* and 
Q = P/mv*. After tedious but straightforward calcula- 
tion, one arrives to the following expression for the mod- 
ulus of the dimensionless relative momentum 



the K-difference introduced in ref. [Tg- Clearly Q(qi, q 2 ) 
represents the modulus of momentum of the first (second) 
particle in the rest frame of second (first) particle. 

We suppose now that the two particles in the frame S 
have momenta qi = q + dq and q 2 = q respectively and 
calculate the value of Q(q+dq, q). One easily obtains 

Q(q + dq,q) = d {K) q , (5.7) 

being q = |q| and d {K) q given by Eq. i|4.7|) . 

The physical meaning of ^-differential d {K} q immedi- 
ately follows. The modulus of the infinitesimal difference 
of the two particle momenta dq, in the frame S, becomes 
d{ K }q, if this difference is observed in the rest frame of one 
of the two particles. The meaning of the K-dcrivativc fol- 
lows also readily. If d/dq is the derivative in the frame S, 
the deformed derivative d/d^q represents an ordinary 
derivative in the rest frame of one of the two particles. 



VI. RELATIVISTIC SUMS AND k-DEFORMED 
LOGARITHM AND EXPONENTIAL 

The ordinary logarithm h(x) = hi(x) is the only ex- 
isting function (unless a multiplicative constant) which 
results to be solution of the function equation h{x\x 2 ) = 
h(x\) + h{x-i)- Let us consider now the generalization 
of this equation within the special relativity obtained by 
substituting the ordinary sum by the K-sum of the di- 
mensionless relativistic momenta 



h(xix 2 ) = h(xi) © h(x 2 ) 



(6.1) 



We proceed by solving this equation, which assumes the 
explicit form 



h(xix 2 ) = h{x\) \j 1 -j- k 2 h(x 2 ) '■ 



+ h{x 2 ) Vl + K 2 ^^ 



(6.2) 



After performing the substitution h(x) = kT x sinh ng(x) 
we obtain that the auxiliary function g(x) obeys the 
equation g{x\x 2 ) = g(x±) + g{x 2 ), and then is given by 
g(x) = Alnx. The unknown function h(x) becomes 



h(x) = 



sinh(Klna:) 



H 



(6.3) 



where we have set A = 1 in order to recover, in the limit 
k — > 0, the classical solution h(x) = ln(x). The function 
given by Eq. 1)6.3(1 in the following is denoted by ln {K} (x) 
and can be written also in the form 



Q(qi,q2) 

being 



\ 



/ K \ 2 2 
( qi e q 2 ) - k 2 (qi x q 2 ) 

1 - k 4 (qi x q 2 ) 2 



qi q 2 = (Wl + K 2 ql - q 2 Jl + K 2 q 2 , (5.6) 



(5.5) 



2 k 



(6.4) 



Holding exp {K} (ln {K} (x)) = ln {K} ( exp {K} (a;)) = x it re- 
sults that ln {K} (a;) is the inverse function of exp {K} (x) 
and thus can be viewed as a generalization of the or- 
dinary logarithm, ln(a:) = ln {0} (x), in the framework of 
special relativity. 
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It is straightforward to obtain the exp {K} (x), starting 
directly from K-sum of the dimensionless relativistic mo- 
menta, as solution of the functional equation 

K 

which in the classical limit reduces to the functional equa- 
tion exp(xi) exp(x2) = exp(xi+a;2) defining the ordinary 
exponential. 

In the Appendix, the main mathematical properties 
of the functions K-exponential and K-logarithm are re- 
ported. 



VII. MAXIMUM ENTROPY PRINCIPLE 

We start by recalling some elements of classical sta- 
tistical mechanics. We consider the Boltzmann-Gibbs- 
Shannon entropy 



s = -jd 3 p /(p)Hf(p)) , 

and the constraints fuctional 



C — a. 



d 3 pf(p)-l 



M-Jd 3 p g(p)/(p) 



(7.1) 



(7.2) 



where the constants a and a = {a t , a 2 , ai} are the 
I + 1 Lagrange multipliers, while the I moments M = 
{Mj, M 2 , Mi} are the mean values of the I functions 
Sip) = {gi(p)i g 2 (p)> ■••) Slip)}- The variational equation 



Sf 



(S + C) = 



(7.3) 



means that the entropy S(p) is maximized under the con- 
straints imposing the conservation of the norm of f(p) 
and the a priori knowledge of the values of the I moments 
Mj of f(p), namely 



d 3 pf(p) = l 



Mi 



d 3 P Sjip)f(p) • (7-4) 



The maximum entropy principle expressed by the varia- 
tional equation (|7.3|) yields the following statistical dis- 
tribution 



Hp) = -exp(-a-gO) 
e 



(7.5) 



In a classical many body particle system the collisional 
invariants are essentially the particle number and the 
kinetic energy W{p) = p /2m. Then we need to in- 
volve two Lagrange multipliers a = /i/fc^T, a I = 
1/fcgT and one moment Mi = J d 3 p gi(p) f(p) — 
J d 3 p W{p) f(p) —<W>. In this way the distribution 
function given by Eq. 1)7. 5|) reduces to the Maxwell- 
Boltzmann one of classical statistical mechanics 



(7.6) 



Let us consider now the following generalized trace 
form entropy 



S 



d 3 P /(p)A(/(p)) 



(7.7) 



being A(/) a function generalizing the ordinary loga- 
rithm. The variational equation l|7.3|l with S and C given 
by H7.7J) and 1)7. 2|) respectively, produces the distribution 
function f — f{p) defined through 



A /A (/) = -a-g(p)+a 



(7.8) 



The maximum entropy principle, in the case of the 
ordinary statistical mechanics, relates the Boltzmann- 
Gibbs-Shanon entropy and the Maxwell-Boltzmann dis- 
tribution by means of a twofold link. Firstly, the distribu- 
tion is obtained by maximizing the entropy under proper 
constraits. Secondly, both the entropy and the distribu- 
tion are expressed in terms of the same function in direct 
and inverse form. Indeed, the Maxwell-Boltzmann distri- 
bution is given in terms of the ordinary exponential while 
the Boltzmann-Gibbs-Shanon entropy is defined starting 
from the ordinary logarithm. Clearly, this is a very strong 
interpretation of the maximum entropy principle. In gen- 
eralized statistical theories, customary, is used a more 
weak interpretation of the maximum entropy principle 
by relaxing the second link. 

Our next task is to use the strong interpretation of 
the maximum entropy principle holding for the ordinary 
statistical mechanics in order to generalize minimally it. 
The requirement of the twofold link between entropy and 
distribution function also in generalized statistical theo- 
ries imposes that the function A(/) is a solution of the 
following functional-differential equation 



!/A(7) 



AA(//a)+r? 



(7.9) 



being {a, A, 77} three real constants. Indeed, the distri- 
bution given by 17. 8f) assumes now the form 



f{p) = ahr 



sip) - 



A 



(7.10) 



being A~ 1 (x) the inverse function of A(x). Note that 
the constants A and i] introduce a scaling in the Lan- 
grange multipliers. Clearly in the classical limit it must 
result: A(x) — > ln(x), A _1 (x) — > exp(x) and {a, A, 77} — > 
{l/e,l,0}. 

The problem related to the determination of the gen- 
eralized logarithm A(/) satisfying Eq. (|7.9|) and the con- 
ditions A(l) = 0, A'(l) = 1 and 0A(0) = 0, can be solved 
easily. Beside the ordinary logarithm, Eq. I|7.9|l admits 
other new solutions. The general solution of this equa- 
tion defines a three-parameter family of generalized log- 
arithms indicated by A(/) = ln {K r f ,(/) being 



hi 



>(/) 



j?r-\-K. K JT — K 



, (7.11) 
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while {re, r, are three free parameters. The constants 
a, A and r\ are given by 



A = 



1 + r — re 



[1 + 7* — K) 



(l + r + k 
rj = (A-l) 



(k + r)<; K + (k - r)s K 



(7.12) 
(7.13) 
(7.14) 



Some specific choice for the parameters {re, r, c}, gener- 
ating particular expressions for the function A(/) and for 
its inverse function A -1 (a;) = exp (x), have been al- 
ready considered in mathematical statistics, information 
theory and statistical mechani cs [2CM. As stressed previ- 
ously in ref.s [II H H EJ 122, l23j . it is shown that 
it is possible to develop selfconsistent statistical theo- 
ries starting from any generalized entropy satisfying some 
standard requirements. Fortunately, all these require- 
ments are satisfied by the entropy defined through (|7.7|) 
withA(/)=ln {Krs} (/). 

Let us focus now our attention to the solution (|7.11() 
corresponding to the particular choice {k, r = 0, tr = 1} 
[l8| . It is easy to verify that results In Q 1} (x) = 
In, * i (x) and exp „ , (x) 



exp {K} (x) with 



2k 



(7.15) 



exp {K} (x) = (^\/ 1 + k 2 x 2 + Kx^j . (7.16) 

Thus the re-exponential and re-logarithm, emerging in the 
special relativity in the place of the ordinary exponen- 
tial and logarithm, reappear also in the statistical theory 
through the maximum entropy principle. It is immediate 
to verify that the reality and positivity of the parameters 
a and A impose that — 1 < k < 1. 

Furthermore, in the case {re, 0, <j}, Eq. (|7.11() produces 
a two parameter generalized logarithm and exponential 
which results to be the scaled re-logarithm and the scaled 
re-exponential respectively, defined as 



(7.17) 



exp 



^x) = ^ w (xy/u^^+]n iK ^) . (7.18) 



We remark that the functions exp (x) and \n {K} (x) 
generalizing the ordinary exponential and logarithm re- 
spectively, emerge both in special relativity and indepen- 
dently in many body physics, through the maximum en- 
tropy principle, the cornerstone of statistical mechanics. 
This interesting fact suggests us to consider with a par- 
ticular attention the statistical theories, enforced by the 



entropies related to the generalized logarithms defined 
through (|7?T3|) and lf7T7jl . 

In the next section we consider the statistical theory 
involving the generalized logarithm (|7.15|l . 



VIII. GENERALIZED DISTRIBUTION 
FUNCTION 



Let us consider the entropy S 



< ln {K} / > related 



to the generalized logarithm defined through H7.15(l . Here 
the distribution function / = f(x,p) depends on the four- 
vectors position x and momentum p. The explicit form 
of this entropy is given by 



d 3 P /ln w (/) 



(8.1) 



and after maximized under the constrains (|7.2II by solving 
the variational equation 



jj(S + C) 







yields the distribution function 



/ = a exp 



{«} 



with g = g(x,p) and 



1 - re 



1/2k 



1 + K v 

A = y/l - re 2 , 



(8.2) 
(8.3) 

(8.4) 
(8.5) 



being |re| < 1. Before to proceed further we explain 
briefly the meaning of the constants a and A, related 
by 1/A = In (1/a). We write the entropy in the form 
S = f n d 3 p a (/) where a (/) = -/ ln {K} (/) is the en- 
tropy density. It is easy to verify that for f — a the en- 
tropy density assumes its maximum value, which is given 
by a m ax = a/ A. 

We focus now our attention on the form of g = g(x,p). 
In ref. p7| it is shown that for a relativistic many 
body system, in presence of external electromagnetic 
fields, the more general microscopic invariant has the 
form {p v + qA u /c) U,y+ constant, being U v the hydro- 
dynamic four- vector velocity with U V U V = c 2 . Then we 
can pose <?i = (p u + qA^/c) U v — mc 2 , a x = 1/ UbT and 
ao = /i/fc^T. The distribution l|8.3|) . in the case of rela- 
tivistic statistical systems, assumes the form 



/ = a exp 



{p v + eA v /c) U v -mc 2 - n 
ALT 



(8.6) 



and results to be quite different with respect the Jut- 
tner distribution of the ordinary relativistic statistical 
mechanics, where in plac e of re-exponential appears the 
ordinary exponential |27j. The distribution 18.6*1 . in the 
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global rest frame where U v — (c, 0,0, 0) and in absence 
of external forces {A v — 0), simplifies as 



f = a exp 



,: \ Xk B T 
being W the relativistic kinetic energy 



•7) 



W = y/m 2 c A + p 2 c 2 - mc 2 , (8.8) 

and in the k — > limit reduces to the classical distribu- 
tion (|7.6|l . It is remarkable that the above distribution 
presents a power law asymptotic behavior, namely 



/ cx W~ 1/K , (W -> oo) , 



(8.9) 



in contrast with the ordinary relativistic distribution 
(Maxwell-Boltzmann distribution where the energy is 
given by its relativistic expression) which decays expo- 
nentially. 



IX. GENERALIZED STATISTICS AND 
KINETICS 

In this section we report the main features of statistical 
theory, based on the entropy JH3J, some of which are 
already known in the literature. 

Equiprobability: The non negative entropy (|8.1I) 
achieves its maximum value at equiprobability, f(p) = 
I/O, for Vp, and this value is S = ln, K ,fi. 

Thermodynamic stability: The entropy (|8.1(l is concave 



S[tfi + (1 - t)h] > tS[h] + (1 - t)S[f 2 ] 



(9.1) 



so that the system, in thermodynamic equilibrium, is sta- 
ble E3- 

Lesche stability: The Lesche stability condition is sat- 
isfied by any physically meaningful quantity. In partic- 
ular it is a necessary condition to experimentally detect 
a physical observable. The entropy given by l|8.1|l . being 
a physical quantity depending on a probability distribu- 
tion, should exhibit a small relative error 



R 



S[f] - S[f] 



su P (S[/D 



(9.2) 



with respect to small changes of the probability distribu- 
tion 



= 11/ -/'II 



(9.3) 



Mathematically this implies that: For any e > there 
exists 5 > such that R < e holds for all distribution 
functions satisfying D < 5. The Lesche stability con- 
dition holds for the Boltzmann-Gibbs-Shannon entropy 
and is proved also for the entropy (|8.1|) in ref. [2Sl. In 
the thermodynamic limit it is obtained 6 — ^/(i-M). 

Relativistic kinetic equation: The statistical distribu- 
tion (|8.7|) can be viewed as stationary case of a distribu- 
tion function / = f{x,p) 1 describing a relativistic many 



body system, governed by the generalized Boltzmann 
evolution equation [Tg| 



p v d w f-mF v 



„ df f d 3 p' d 3 Pl d 3 p[ 



dp- J p'° p° pf 
o(/'®/D-o(/®/i) 



G 



(9.4) 



In the latter equation the standard notations of the ordi- 
nary Boltzmann kinetics are used and a(f) is a positive 
and increasing arbitrary function, which does not affect 
the stationary form of /. The factor G is the transition 
rate, depending only on the nature of the two body parti- 
cle interaction. The composition law defined through 



(9.5) 



is a one-parameter generalization of the ordinary product 
having the following properties 

i) associative law: (/ h) ® w = f ® (h <8> w), 

ii) neutral clement: / <g)l = l<g)/ = / ! 

K K 

iii) inverse element: / ® (1//) = (1//) <8> / = 1, 

iv) commutative law: / ® h — h <g> / '. 

K K 

Furthermore it results / (8>0 = 0(8)/ = and the general- 

K K 

ized division can be defined through / 0/i = /<g>(l/ft). 

H-theorem: The four-vector entropy S u — (S°,S) it 
was defined as 



(9.6) 



The temporal component S° of this four-vector coincides 
with the K-entropy defined previously through Eq. i|8.1|) 
while the spatial component S is the entropy flow. Start- 
ing from the evolution equation l|9.4l) it has been demon- 
strated the H-theorem, which represents the second law 
of thermodynamics and states that the entropy produc- 
tion is never negative, and, in equilibrium conditions, 
there is no entropy production. Namely, the following 
relation was obtained 



d v S" > 



(9.7) 



which is the local formulation of the relativistic H- 
theorem and represents the second law of the thermo- 
dynamics [l9l |. 

Entropy transformation: After recalling the identity 
d 3 p/p° = d 4 p 26(p°) Sip^Pf, - m 2 c 2 ) the definition i|9~B|) 
assumes the form 



d 4 p 2 6(p») *(p% - mV) p v f ln {K} / 



(9.8) 

In this expression d 4 p is a scalar because the Jacobian of 
the Lorentz transformation is equal to unit. Then since 
p v transforms as a four-vector, we can conclude that S v 
transforms as a four-vector like the particle four-flow 



N u = d 4 p 2 0(p°) Sfjfpn - m 2 c 2 ) p v f 



(9.9) 
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for which the conservation law of the total particle num- 
ber imposes 



d v N v = 



(9.10) 



Calculation of the parameter k: We focus now atten- 
tion to the general discussion of sec. II. The results ob- 
tained in sect. II permit us to consider the problem re- 
lated to the calculation of the parameter k from a very 
different and more physically sound point of view with 
respect the one of ref. [l!j. Clearly, the velocity in- 
troduced through | At | = v*/c does not emerge within the 
one-particle relativistic theory. In general v* could de- 
pend on c namely v* = v*(c) with the condition that 
f*(oo) < oo. The fact that the condition \k\ < 1 and 
then i>* < c, follows from the maximum entropy prin- 
ciple, supports the supposition that u* could emerge in 
the framework of the many body theory. Taking in mind 
that in statistical mechanics the only emerging velocity is 
the thermal velocity, in the following wc identify z;* with 
the thermal velocity of the many body system. 

The particle thermal energy, being the mean particle 
kinetic energy (W), is related to the thermal velocity 
u* = kc according to the relativistic relationship 



(W) = 



1 



i 



(9.11) 



Clearly, in a statistical theory, the mean particle kinetic 
energy (W) is considered a known quantity. Then the 
expression of k in terms of (W) follows immediately 



K = 1 



(9.12) 



so that the theory does not contain free parameters and in 
the classical limit c — > oo, results k = 0. From i|9.12[l fol- 
lows that lim^oo kc = ^2 < W > jm < oo as required 
from the theory developed in sect. II. Note that in the 
case we know experimentally the distribution function, 
the value of k and consequently, by using (|9. 12|) , the value 
of the ratio (W) /mc 2 can be obtained from the slope of 
the distribution tails in a log-log plot. The normalization 
condition for the distribution function / and the defini- 
tion of the particle mean kinetic energy, starting from /, 
permit us to determine univocally T and /i . In this way 
one can link (W) and then n to the temperature of the 
system. 

Bearing in mind the meaning of the parameter k we 
rewrite the distribution function 1|8.7[) in the form 



/ = a exp 



W -n 
k a T' 



(9.13) 



being W the relativistic kinetic energy while the temper- 
ature T' is given by 



r 



T 



(9.14) 



We remark that the temperature T' is contracted with 
respect to the temperature T by a factor which, curiously, 
has the same form of the one entering in the definition 
of the relativistic temperature suggested by Einstein p9| 
and Planck |30j . 



X. APPLICATIONS 

In the last four years the present theory has been used 
successfully to study some systems which manifestly can- 
not be treated within the ordinary statistical mechanics. 
Clearly, the most natural applications regard relativistic 
many body systems. 

The first application, considered in ref. [19|, concerns 
the cosmic rays. Since long time it is known that the 
cosmic rays spectrum, which extends over 13 decades in 
energy, from a few hundred of MeV (10 8 eV) to a few 
hundred of EeV (10 20 eV) and spans 33 decades in parti- 
cle flux, from 10 4 to 10 -29 [m 2 sr s GeV] -1 , is not expo- 
nential and, then, it violates the Boltzmann equilibrium 
statistical distribution cx exp(-E/k B T) [Hill El- 0n 
the other hand, it is known that the particles compos- 
ing the cosmic rays are essentially the normal nuclei as 
in the standard cosmic abundance of matter. Then, the 
cosmic rays can be viewed as an equivalent statistical 
system of identical relativistic particles with masses near 
the mass of the proton (938 MeV). These characteristics 
(relativistic particles with a very large extension both for 
their flux and energy) yield the cosmic rays spectrum an 
ideal physical system for a preliminary test of the correct- 
ness and predictability of any relativistic theory. In ref. 
|l9j , a high quality agreement between the predictions of 
the present theory and the observed data, in the whole 
cosmic ray spectrum, has been found. This agreement 
over so many decades is quite remarkable. 

Other physical applications of the theory regard the 
formation of quark-gluon plasma |3l| and kinetic models 
describing a gas of interacting atoms and photons |32j . 

On the other hand, the theory has been applied suc- 
cessfully also to the study of natural or artificial systems 
exhibiting a limiting velocity in the propagation of the in- 
formation (like the one imposed by light speed in physics) 
where a mechanism analogous to the relativistic one can 
emerge, deforming the distribution function. In ref. |33| 
the problem of the fracture propagation has been con- 
sidered within the present theory. Others applications 
have been considered in games theory [34| . in economy 
for study the income distribution |35j and in construct 
financial models [3(j, etc. 

We come back to the high energy physics being the 
more proper field of application of the theory. In the 
following, as further application, we consider the open 
problem of the black hole physics regarding the intrigu- 
ing question of the Bekenstein-Hawking area law p^.ll5|. 
This law asserts that the entropy of a black hole scales as 
the area A of the event horizon namely S oc A. It is well 
known that the ordinary Bolzmann statistical mechanics 
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and thermodynamics fails to justify this law and in con- 
trary predicts that the entropy scales as the volume V 
of spatial region delimited by the event horizon, namely 
SocV. 

Indeed, if we indicate with N the degrees of freedom or 
microscopic components in the volume V and with M the 
number of states accessible by each component, within 
the Boltzmann statistics the number of microstates of 
the whole system is given by fi oc M N with N >> 1. 
After introducing the notations 

7i(0)=cxpiV , m(0) = mM , (10.1) 

the number of microstates can be written also in the form 

fi oc n(0) m(0) . (10.2) 

Customary one assumes a uniform distribution of the 
components in the volume V obtaining N oc V. Thus 
the entropy 



S = lnSl , 
of the black hole results to be given by 
SocV . 



(10.3) 



(10.4) 



This result of the Boltzmann statistics, already known 
in the literature, is manifestly in contradiction with the 
Bekenstein-Hawking area law. 

Very recently, in ref. [37j, has been suggested to treat 
this problem within non-standard statistics. Even if the 
authors limit their discussion to the nonextensive |3£j | as- 
pects of the problem the proposed approach has a general 
validity and can be adopted also in the present case. 

Clearly the black hole is a relativistic system and ap- 
pears natural to write its entropy in the microcanonical 
picture as 



5 = ln w n 



(10.5) 



Furthermore the present formalism generalizes Eqs. 
(|TTH|) and ltTIO|l as follows 

n(«) -exp {K} iV , m{ K )=ln {Kj M , (10.6) 



We remark that the power law behavior of the various 
physical quantities is a feature of the present theory in 
contrast with the exponential and logarithmic behaviors 
in classical statistical mechanics. Finally, we assume, as 
customary, a uniform distribution of the components in 
the volume V, namely N cx V and obtain the following 
expression of entropy 



After posing 



S cx V m{K) 



m(n) = 2/3 



(10.12) 



(10.13) 



the Bekenstein-Hawking law, namely S oc A, follows im- 
mediately and thus we can conclude that the present sta- 
tistical theory is consistent with the area law of black hole 
entropy. This result is in agreement with ref. |3?1 | where 
the authors propose a relationship analogous to the one 
given by Eq. (|10.13|) . 



XI. EPIMYTHION 

We have shown that within the standard framework of 
the special relativity, the Lorentz transformations lead 
to a proper one-parameter generalization of the entropy 
just as for other physically meaningful quantities (e.g. 
momentum, energy, etc). 

The new entropy generates a coherent and selfconsis- 
tent relativistic statistical theory (both statistical me- 
chanics as well as kinetic theory) preserving the main fea- 
tures (maximum entropy principle, thermodynamic sta- 
bility, continuity condition, H-theorem, etc) of the ordi- 
nary statistical theory which riemerges in the classical 
limit. 

The distribution function predicted by the present the- 
ory is a one-parameter, continuous deformation of the 
classical Maxwcll-Boltzmann distribution exhibiting a 
power law tail in accordance with the experimental evi- 
dence in several relativistic many body systems. 



Q oc n( K ) m(K) 



(10.7) 



After recalling the property (|12.15l) the number of mi- 
crostates becomes 



ft cx exp , , s (m(K)N) , 

{K/m(K)} v v ' ' ' 

while the black hole entropy can be written as 



S cx In 



exp (m(K)N) 



(10.8) 



(10.9) 



In the thermodynamic limit N >> 1, then Eqs. 
and (|10.9|) predict that the number of microstates and 
the entropy scales as power laws 



n cx iv m(K)/K 

ScxN m ^ . 



(10.10) 
(10.11) 



XII. APPENDIX 

We start by considering the K-sum for the dimension- 
less relativistic momenta 



and after posing t\ — t + dt and t% — —t we can introduce 
K-differential through 



d {K} t =(t + dt)et 

= (t + dt) e (-t) 

dt 



VI + nH 2 



(12.2) 
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The K-derivative d/d {K} t admits as eigenfunction the k- 
exponential, namely 



exp {K} (t) = exp {K} (t) , (12.3) 



being 



exp {K} (t) = (Vl + « 2 * 2 + «t) 1A , (12.4) 
exp fK ,(i) = exp f — arcsinh(zcf) J . (12.5) 



a typing error (the first minus in the denominator of the 
right hand side must be replaced by a plus so that after 
some simple algebra Eq. (5.31) of |[l.9] reproduces the 
present result). 

The K-logarithm ln {K} (f), defined as the inverse func- 
tion of the of ^-exponential, namely ln {K} (exp t) = 

ex p {K) ( ln {K> *) = *> is s iven b y 



t K - t~ K 
2k 



ln {K} (i) = — sinh(/iln£) 



(12.21) 
(12.22) 



The K-exponential can be viewed as a generalization of 
the ordinary exponential 



and can be viewed as a generalization of the ordinary 
logarithm 



exp {0> W = exp W ' 



and has the properties 



(12.6) 
(12.7) 



%)>(*)=M*) » 



(12.23) 
(12.24) 



exp w (t)GC~(R) , 


(12.8) 


ln {K} (i)eC 00 (R+) , 


^ex P{K} (i)>0 , 


(12.9) 


Al» w (*)>0 , 


d 2 


(12.10) 




^2 ex P w W>0 , 




ex P {K} (-°°) = 0+ , 


(12.11) 


ln w (0+) = -oo , 


exp w (0) = l , 


(12.12) 


V}( 1 ) = ° ' 


exp {K} (+oo) = +oo , 


(12.13) 


ln { «}(+ 00 ) = +°° > 


ex P{K} (t)ex P{M H) = l . 


(12.14) 


ln w (l/t) = -ln w (t) 



Furthermore it results 

( ex P { «>(*)) r = ex P {KM ( rt ) ' ( 12 - 15 ) 
ex P{K} (ti) ex P{K} (ta) = ex P{K} (h © t 2 ) , (12.16) 
exp {M (t x )® exp {K} (*a) = exp {K} (t x + t 2 ) ,(12.17) 

where the product <g> is defined through Eq. 119.51) . The 
following asymptotic relations also hold 



The K-logarithm, just as the ordinary logarithm, has the 
properties 

(12.25) 
(12.26) 

(12.27) 

(12.28) 
(12.29) 
(12.30) 
(12.31) 

Moreover it has also the properties 

m w (r)=rln {rK} (t) , (12.32) 

ln {K} (ti i 2 ) = ln {K} (ti) © ln {K} (i 2 ) , (12.33) 
ln {K} (ti® ta) = m w (ti) + ln {K} (i 2 ) , (12.34) 

while the following asymptotic formulas hold 

K 2 \ t 3 



exp^-l + t+^ + il-n 2 )^ , (12.18) 



exp (t) ~ ^/stl^/l"! . (12.19) 

{M t->±oo 

Finally, after tedious but straightforward calculations the 
following formula can be obtained 



,r-l , ,x |2<c|~ r r fe _ V 

r exp (—t)dt= — — -, 



T(r) 



(12.20) 

holding for r|«| < 1. Note that this integral is given 
also in [l9| by Eq. (5.31) which unfortunately contains 



ln - (1+t) t ^4 + ( 1 + y) 

Furthermore holds the very useful formula 
and finally one obtains [l^ 



di In 



\2k\ 



t J 1 + (r 



r— 1 



(12.35) 

(12.36) 
(12.37) 

(12.38) 

- F(r) . 
(12.39) 
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It results evident that the above used procedure, to in- 
troduce the ^-exponential and ^logarithm, starting from 
the additivity law of the momenta, can be used to obtain 
new functions if we start from other generalized additiv- 
ity laws. For instance, starting from the additivity law 
of the velocities we have d {K} t = (1 — K 2 i 2 ) _1 cft, while 
instead of ^-exponential we obtain the function 

^>=(^0" 2 "' (12 - 40) 

= exp ^-arctanh(frf)^ , (12.41) 

whose inverse function assumes the form 

1 t K — t~ K 
V W =- 7^-7-, > ( 12 - 42 ) 

K t K +t K 

= -tanh(Klnt) . (12.43) 

K 

The function ip(t) G C°°(I) with I =]— |/c| _1 , 



is connected to exp {K) (t) through 

" (t) = exp »(7r^F) ' {l2M) 

At this point one could be tempted to replace exp {K} (i) 
with the function (p(t) as the right one for generalize the 
ordinary exponential. 

There are several reasons to consider exp {K} (t) as the 
more proper generalization of the ordinary exponential. 
Beside the motivations emerging within the one-particle 
relativistic dynamics, there is the one related to the max- 
imum entropy principle. This principle, the cornerstone 
of statistical physics, selects without any ambiguity (see 
sect. VII) just the K-exponential as the more natural gen- 
eralization of the ordinary exponential. More precisely 
the function exp {K} (t) is the only one existing which 
emerges simultaneously both in the one-particle relativis- 
tic dynamics as well as in the many-body physics through 
the maximum entropy principle. 
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